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ABSTRACT 

Context. Monitoring of the photometric and chromospheric HK emission data series of stars similar to the Sun in age and average 
activity level showed that there is an empirical correlation between the average stellar chromospheric activity level and the photometric 
variability. In general, more active stars show larger photometric variability. Interestingly, the measurements and reconstructions of 
the solar irradiance show that the Sun is significantly less variable than indicated by the empirical relationship. 
Aims. We aim to identify possible reasons for the Sun to be currently outside of this relationship. 

Methods. We employed different scenarios of solar HK emission and irradiance variability and compared them with available time 
series of Sun-like stars. 

Results. We show that the position of the Sun on the diagram of photometric variability versus chromospheric activity changes 
with time. The present solar position is different from its temporal mean position as the satellite era of continuous solar irradiance 
measurements has accidentally coincided with a period of unusually high and stable solar activity. Our analysis suggests that although 
present solar variability is significantly smaller than indicated by the stellar data, the temporal mean solar variability might be in 
agreement with the stellar data. We propose that the continuation of the photometric program and its expansion to a larger stellar 
sample will ultimately allow us to constrain the historical solar variability. 
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1. Introduction 

The comparison of the Sun with Sun-like stars helps to improve 
our understanding of both solar and stellar physics. One of its 
important applications is the study of the solar irradiance vari- 
ability. The consecutive measurements of the solar irradiance 
have only been available for the last three decades. The pho- 
tometric observations of the large number of Sun-like stars can 
help to extend the time coverage and constrain the historical vari- 
ability of the solar irradiance. 

The first monitoring of the photometric flux and chromo- 
spheric HK emission (which characterises the amount of non- 
thermal heating in the chro mosphere and consequently activ- 
ity, see iBaliunas et al.| [T995) of the stars similar to the Sun in 
age and average activity level revealed an empirical correla- 
tion between the stellar photometric variability and mean chro- 
mospheric activity level (Lockwoodet al. 1992). The continua- 
tion of the program (see e.g-J Lockwood et al.ll997t iRadick et al.l 
ll998tlLockwood et al.l2007l) showed that the variability of a star 
with solar chromospheric activity level is expected to be about 
0.1 % in the Stromgren filters b and y (centred at 467 and 547 
nm, respectively). At the same time the regression of the ob- 
served total solar irradiance (TSI) variability to the variability in 
Stromgren filt ers b and y yielded a somewhat smaller value of 
about 0.04 % (IRadick et alJll998l) . 

Several a ttempts were made to explain this disagreement. 
For example, iFoukall (1 19941) proposed that the relatively low so- 
lar variability can be explained by the dependenc e of the f aculae 
to starspot ratio on the magnetic activity. Later, Radick (1994) 
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show ed that this ex planation is not consistent with the stellar 
data. Sch attenl d 1993b suggested that the special position of the 
Earth-based observer who sees the Sun from its equatorial plane 
can reduce the variability of solar brightness. However, detailed 
calculations showed that the equatorial position of the observer 
leads to only a 30% (most probable value) decrease in the solar 
variability and does not affect the m ean level of the chromo- 
spheric activity dKnaack et al.ll200H) . The out-of-ecliptic mea- 
surements of the Sun are needed to validate the model and ac- 
curately assess t he effects from the hemispherical asymmetry of 
the density flux dVieira et al.ll2012l) . 

A significant part of the disagreement can also be at- 
tributed to the s election effect in the stellar sample. For example, 
Hall et a using another stellar sample, presented a dif- 

ferent version of variability versus activity regression. It is based 
on the monitoring of 28 Sun-like stars over time period s rang- 
ing from 6 to 1 5 yr. The slopes of the lLockwood et aTl {2007) 
and lHall et alJ (2009) regressions are approximately the same, 
but the lHall et alJ (12009) regression indicates about 0.1 dex less 
variability. 

Recently Judge et all (1201 2l) proposed comparing the recon- 
structed time series of the solar irradiance with the photometry 
of Sun-like stars. In this paper we follow up on this idea. We re- 
evaluate the value of the solar variabi lity in Stromgren b and y 
filters obtained bv lRadick et al.l (1 19981) and show that the present 
location of the Sun on the variability versus activity diagram is 
anomalous. 

The continuous measurements of solar irradiance have been 
available only since the beginning of satellite observations (i.e. 
since 1978). Meanwhile, the proxy data indicate that the Sun 
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Fig. 1. Photometric variability versus chromospheric activity. 
The blue stars indicate the stars with observed variability from 
lLockwood et al.l d2007l) . Magenta triangles indicate lower lim- 
its for stars with unconfirmed variability set at about 0.5 be- 
low the regression line. The solid (5.892+1.8101og R' HK ) and 
dashed (7. 180+2. 108 log R^ K ) lines are regressions calculated 
excluding and including stars with unconfirmed variability, re- 
spectively. The shaded area corresponds to the one-sigma er- 
ror. The four red stars indicate HP 10307 HD951 28, HD168009, 
and HD146233 (18 Sco) from lHall et all (l2009h . Drop lines in- 
dicate the correction for the variability of the comparison stars. 
The orange symbols denote the position of the S un according to 
moder n reconstructions (solar symbol) and the lLockwood et alJ 
(2007) estimate (square symbol). 



has been at a maximum p lateau- state in its long- t erm evolution 
durin g the last 50 years (lLockwood et al.l fT999t ISolanki et al. 
2004). Therefore, the amplitude of the 11-year activity cycle is 
expected to be relatively large, while the amplitude of the sec- 
ular compone nt, whose existen ce is actively speculated in the 
literature (cf. iLockwoodl 120 1 ll) . is very small in the most re- 
cent past. At the same time cosmogenic isotope records indi- 
cate that during the last millennia the level of the long-term 
solar activity has been sign ificantly variable dMcCracken et alJ 
2004; Vonmoos et al.1 12006). This implies that the amplitude of 
the 1 1 -year cycle also has been changing in time and, addition- 
ally, that the Sun undergoes periods when the long-term photo- 
metric variations could be relatively strong. Therefore, the tem- 
poral mean of solar variability might be significantly different 
from the presently measured value. To investigate this possibil- 
ity we locate the Sun on the variability versus chromospheric ac- 
tivity diagram, employing different scenarios of solar long-term 
variability. 



2. Comparison of the stellar and present solar 
variability 

Figure Q] presents the regression of the stellar variabil- 
ity versus chromosphe ric activity calculated with data from 
lLockwood et al.l (|2007). To calculate the regression line we 
only consider t he st ars which were deemed variable by 
lLockwood etall d2007h . 

One can see that even for the stars with the same chromo- 
spheric activity level there is some divergence in the photomet- 
ric variabilities. The variability of the less variable stars with 



low magnetic activity would be very difficult to detect (e.g. two 
of the five stars with unconfirmed variability had log R^ K < -5). 
Therefore, the regression was made using only the stars with 

1o 8 R hk > - 5 - 

Orig inally the Lowell O bservatory program contained 41 
stars (lLockwood etal.1 119971) . Some of the stars with uncon- 
firmed variability were later excluded from the program so the 
final set contained only 32 stars. Seven stars from the original 
dataset had log R^ K > -5 but unconfirmed variability. The in- 
ability to detect the variability of the star is either because its 
variability is indeed small or because the comparison stars are 
too variable. In the second case, the exclusion of the star is justi- 
fied; however, in the former case the exclusion of such star from 
the analysis would affect the variability versus chromospheric 
activity dependency. To estimate the possible effect of such stars 
on the regression, we located the stars (magenta triangles in 
Fig. Q]) with unconfirmed variability below the regression line 
on the distance equal to the maximum deviation of the stars with 
confirmed variability from the regression line. While the original 
regression indicates that the root mean square (RMS) variability 
of a star with the solar chromospheric activity level should be 
about 0.001 mag, the new regression yields a value of 0.00072 
mag. These two numbers can be considered as lower and upper 
limits of the variability of a star with solar chromospheric activ- 
ity level, yielded from the present stellar s ample. 

To take into account the new data of lHall et al.l d20 09 ) we 
also calculated two regressions: one using the lHall et alJ (I200S 
data corrected for the variab il ity of the comparison stars and 
another using th e Hall et al] d2009l) data combined with the 
Lockwood et all d2007h data. All these different versions of the 
regression yield solar variability within the error bars, esti- 
mated above (0.00072-0.001 mag). In Fig. Q] we show the vari- 
abilities of four stars fro m the ELO DIE top ten sola r analogs 
dSoubiran & Triaudll2004l) taken from lHall et ail (120091) . 

Establishing the present solar position in Fig. Q] is not a 
straightforward task as there are no long-term measurements 
of the spec tral solar variability in the band filters. Radic k~et al.l 
( 1998) and lLockwood et all ((2007) assumed the total solar irra- 
diance and photometric variabilities to be connected as if they 
were caused by a change in the solar effective temperature. 
This resulted in 0.00044 mag estimation of the RMS of annual 
Stromgren (b+y)/2 values. 

The modern physically-based reconstructions of the solar 
irradiance that attribute the solar variability on the 11 -year 
time scale to the compet ition between the dark sunspots and 
bright active region s (e.g. Lean et alj|2005t iKrivova et al.l l2010; 
Shap iro et al.l201 ll) . give smaller values of the present solar vari- 
ability for the last 50 years: 0.00027, 0.00017, and 0.00017 mag 
respectively. We used the mean of these three numbers to locate 
the Sun in Fig.Q] 

Let us note that if instead of using the annual averages of 
the solar irradiance we emulate the stellar observations and con- 
sider the seasonal averages then the calculated solar variability 
can be increased by the rotational cycle. We cal culated the real 
seasonal means using the lLockwood et al1 (l2007) time series and 
estimated that the increase is generally smaller than 40%. 

The measurements obtained by spectral irradiance monitor 
(SIM) on-board of the solar radiation and climate experiment 
(SORCE) satellite show yet another picture of solar variability 
over the last seven years (lHarder et al.ll2009l) . The RMS of the 
annual (b+y)/2 values for 2004 - 2011 period is 0.00029 mag. 
The extrapolation of the SIM annual (b+y)/2 values to the entire 
solar cycle using the annual v alues of the sun spot number results 
in 0.0004 mag for the RMS. iFontenlal d2011l) suggested that a 
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Fig. 2. Upper panel: Annual (b+y)/2 values versus time ac- 
cording to the spectral s olar irradiance (SSI ) reconst ructions to 
the pa st as calculat ed by iKrivova et al.1 d2010l) (red), iLean et aTl 
(120051) (blue) and lShapiro et al.1 (1201 ll) (magenta: thick and thin 
curves are original data and data with removed secular changes, 
respectively). Stars indicates the SIM/SORCE measurements. 
Lower panel: Calculated for the 20 year time intervals RMS20 
values versus time. 

physics-based extrapolation of the SIM data could put the Sun 
even closer to the regression line from Fig. [1] 

3. Comparison of the stellar and historical solar 
variability 

Our main proposition is that the location of the Sun in Fig. Q]is 
not fixed in time. Therefore, instead of using present solar vari- 
ability for the comparison with Sun-like stars, one should use the 
temporal mean of the variability over a time interval long enough 
to reveal the entire range of solar variability. Because of the high 
solar activity the amplitude of the 11 -year cycle has been rela- 
tively large during the last 50 years. So, if the 11 -year cycle is 
the only contributor to the solar irradiance variability then the 
temporal mean of solar photometric variability is lower than its 
present value (see calculations below). This would imply that the 
disagreement between the stellar and solar photometric variabil- 
ities is even more pronounced than presently thought. Instead, 
a secular component in the solar irradiance can significantly in- 
crease the temporal mean of the solar photometric variability. 

The present consensus is that secular changes in solar ir- 
radiance are proportional to secular changes in solar activity 
proxi es. However, there is no agreement on the scaling coeffi- 
cient (|Lean et a l. 2005; Tappi ng et alJl2007llKrivova et al.l feOlO: 
Schri iver et al 1 1201 ll IVieira et all 1201 ll IShapiro et all I20TTI 
Thuil lier et al.l l2012t Hudge et all 12012V To illustrate the effect 
of the secular variations on the solar photometric variability in 



Fig. [2] we plot variations of annual (b+y)/2 values and magni- 
tudes of RMS20 (the RMS20 value for the year X is the RMS 
of (b+y)/2 values in the [X-9 yr,X+10 yr] dataset) for three re- 
cent spectral solar irradiance (SSI) reconstructions to the past. 
Most of the stars used in this study had been observed for ap- 
proximately 20 years, so the RMS20 value can be used for direct 
comparison of the Sun and stars. According to all three recon- 
structions the secular changes in the solar irradiance have been 
very small during the last 50 years. Therefore, during this period 
the solar variability (i.e. RMS20 value) is roughly the same in 
all three reconstructions. At the same time, there were periods 
before 1950 when the long-term solar activity chang ed signif- 
icantly dMcCracken et al J 120041 ISteinhilber et al.ll2008l) and the 
reconstructions with large and small secular components are dra- 
matically different. 

In order not to be limited by any particular reconstruction 
we introduce two parameters, Vu and Vlt, which describe the 
1 1-year and long-term variability, respectively. We define V\ \ as 
RMS variability of solar annual Stromgren (b+y)/2 fluxes dur- 
ing cycles 21 and 22 (the 1976-1996 period) and Vlt as rela- 
tive change of Stromgren (b+y)/2 flux between 22/23 and the 
Maunder minima. In Appendix [A] we describe a simple model 
which allows us to calculate the solar variability over the millen- 
nia as a function of V u and Vlt parameters, u s ing the cosmo- 
genic isotope records dMcCr acken et al] 120041: IVonmoos et al.l 
2006). We show that while our approach does rely on the as- 
sumptions, they do not constrain its applicability, and by choos- 
ing corresponding Vlt and V\ \ parameters we can roughly repro- 
duce very different scenarios of the solar variability published in 
the recent literature. 

Let us note that V\ \ and Vlt values cannot be directly com- 
pared because Vu corresponds to the RMS variation during a 
20-year interval, whereas Vlt describes an absolute change over 
300 years. 

An example for a chosen pair of parameters is illustrated in 
Fig. IB. 21 (available as Online Material) for the time period from 
400 BC to present. 

The mean level of the solar chromospheric activity can also 
change with time. Currently there is no reas on to suggest that 
it underwent significant changes i n the past (Hall & Loc kwoodl 
120041: Uudge & Saaill2007h . ISaari |2006) showed that the mini- 
mum chromospheric activity of stars with solar metallicity cor- 
responds to log R' = -5 .08. We will now consider two extreme 
cases. The first one corresponds to the absence of a secular com- 
ponent in the mean solar chromospheric activity (which does not 
imply the absence of the secular component in the photometric 
variability). So changes of the solar activity index occur only 
because of the evolution of the solar activity cycle. The second 
case (hereafter, case of strong secular changes) corresponds to 
the assumption that the mean solar chromospheric activity can 
reach a minimum boundary value log R' HK = -5.08 during the 
periods of zero modulation potential. In Appendix lAl we show 
how the mean solar chromospheric activity can be reconstructed 
to the past in both of these cases. 

Figure [3] shows the trajectories of the Sun over last 9000 
years in the photometric variability versus chromospheric activ- 
ity diagram for a few different scenarios. The upper-left panel 
represents the scenario without secular variability. The Vlt = 
0.044% (uppe r-right panel) corresp onds to the secular variabil- 
ity reported bv lKrivova et al.l (120101) . During the last thirty years 
the Sun was unusually active so its mean chromospheric activity 
over the last 9000 years is shifted towards lower values relative 
to the present position. The mean value of log R' HK equals -4.92 
without the long-term trend in the chromospheric activity, while 
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Fig. 3. The same as Fig.fT] but with trajectories of the Sun over 
the last 9000 years (orange dots) for different sets of the long- 
term variability Vlt and the 11 -year variability Vn. The red 
square and cross symbols indicate the positions of the present 
and "mean" Sun. 



with strong long-term trend it equals -4.96. The amplitude of 
the 11-year cycle in Stromgren (b+y)/2 averaged over the last 
9000 years is approximately 40% smaller than its present ampli- 
tude, so if the photometric variability is not affected by a secular 
component then the "mean" Sun is located even lower than it 
is now. However, a secular component in the photometric vari- 
ability shifts the "mean" Sun towards larger variabilities and in 
combination with a displacement towards lower chromospheric 
activities moves it to the regression line. 

Let us note that our basic assumption implies that the 
variability of Sun-like stars is also time-dependent (see 
Appendix [B). In combination with the inclination effect and 
measurement errors this can explain the scatter in the activity 
versus variability diagram, which most probably, is real. 

4. Discussion and conclusions 

We have compared the available photometric records of Sun-like 
stars with a parameterized reconstruction of the solar irradiance 
over the last 9000 years. We assumed different scenarios of the 
solar secular variability and showed that if the solar irradiance is 
not affected by secular variations then there is a significant dis- 
agreement between the stellar variabilities and temporary mean 
variability of the Sun. Part of this disagreem ent can be expl ained 
by the selection effect in the stellar sample dHall et al.l2009l) . rel- 



atively small number of stars, and short time span of the stellar 
observations. 

At the same time, we show that the presence of a secu- 
lar component in the solar irradiance can increase the temporal 
mean of solar variability. By assuming that the temporal mean 
Sun obeys the activity versus regression correlation, and employ- 
ing the currently available datasets to constrain the V\ \ parame- 
ter, we try to estimate the magnitude of the secular component. 
The calculations are presented in Appendix [B] While the cur- 
rently available stellar data yield relatively large amplitude of 
the secular component, careful estimation of the involved un- 
certainties shows that a 0.17 % decrease of the solar flux in the 
Stromgren (b+y)/2 filters during the Maunder minimum com- 
pared to the present value is within 95% confidence level. The 
trajectories of the Sun in the photometric variability versus chro- 
mospheric activity diagram which correspond to this amplitude 
of the secular component are shown in the middle panels of 
Fig.|3](for two scenarios of the chromospheric activity variabil- 
ity). We note that the positions of the Sun-like stars on the vari- 
ability versus activity diagram are also not necessarily fixed in 
time, and this can affect the variability versus activity regression. 
This is taken into account in Appendix [B] and therefore, while 
the Vlt = 0.17% value is consistent with stellar data, the cor- 
responding temporal mean position of the Sun is still below the 
regression line. The Vlt = 0.17% value corresponds to 1.9-2.7 
W/m 2 TSI change (see AppendixIBl. 

We note that this number is based solely on the empirical 
relationship between the stellar variability and activity and does 
not depend on mechanisms of stellar and solar variability. It is 
given by the present stellar sample and should be carefully re- 
evaluated when more stellar data become available and the low 
end of the activity sequence is studied. The continuation of the 
stellar photometric program will thus ultimately allow one to 
prove or disprove the existence of the secular variations of the 
solar irradiance. 
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Appendix A: Model description 

The annual solar Stromgren (b+y)/2 flux values F(t) can be de- 
composed into the long-term F hT ( t) and cyclic F n (t) compo- 
nents (see, e.g. [japping et alj|2007h 



F(f) = F LT (0 + F n (t). 



(A.l) 



While cyclic component F n (t) varies on the 1 1 -year time scale 
and controls the activity cycle, the long-term component F hT (t) 
varies on secular time scales and leads to a changeable level of 
the flux during the solar acti vity minima. 

Most of the stars in the lLockwood et all (|2007) dataset are 
observed over an approximately 20-year period. Therefore, to 
compare stellar and solar variabilities we will introduce the 
RMS20 value. For any given year f, it is the RMS variability of 
the dataset which consists of twenty annual flux values (from 
year t-9, till year t+10). 

The total variability of the solar flux can be calculated as the 
geometrical sum of the long-term and cyclic variabilities 

RMS 20 CF(0) 2 = RMS 20 ( J F LT (f)) 2 + RMS 20 CF n (f)) 2 . (A.2) 

We introduce parameter Vu as RMS variability of the an- 
nual (b+y)/2 flu x values for the 1976-1996 period (solar c ycles 
21 and 22) ThelKrivova et all (l2010h . lShapiro et al l (1201 ll) . and 
iLean et al I (12001 " reconstructions show the value of V\ \ param- 
eter equals 0.00019 mag, 0.00018 mag, and 0.000032 mag, re- 
spectively. 

To calculate the cyclic component backwards in time, we as- 
sume that its amplitude is proportional to the long-term solar 
activity (pit). The latter can be reconstructed over the millennia 
from the cosmogenic isotope data, whi ch gives a 22-year run- 
ning mean of the m odulation potential (Steinhi lber et al.l l2008; 
iHerbst et afll2010t> . Thus, the cyclic RMS variability for any 
given year t can be calculated as 

RMS 20 CF n (f)) = Vn-* (0 



<P(to) 



(A3) 



where 0(f) is the modulation potential and <f>(to) is the mean 
value of the modulation potential for the 1976-1996 period. 
The modulation potential is deduced from th e composite of 
the neutron monitor da t a (lUsoskin et al.1 20051) and 10 Be data 
dMcCracken et al. 2004t Vonmoos et al.ll2006l) . It is normalised 
to the lCastagnoli*& LaJ(l980) local in terstellar spectra. The ho- 
mogenisation procedu re is discussed in Shapiro e t al.l(l201 ll) and 
Thuill ier et all d2012l) . The average uncertainty of the mo dula- 
tion potential changes is ~80 MeV dVonmoos e t al. 2006]) and 
does not have a significant contribution to the temporary mean 
of the solar variability. 

To calculate the RMS 2 o(,F LT (f)) contribution we introduce 
parameter Vlt as a relative change of F hT between the present 
and the Maunder minimum. The changes of the solar flux be- 
tween three recent solar minima are within the uncertainties of 
the reconstruc tions and measure ments of the solar irradiance 
(see review bv lErmolli et al.l2012l) . so the F LT component can be 
approximated by a constant for the last solar cycles. Therefore, 
for simplicity we define Vlt as the relative change of the (b+y)/2 
flux between the 1996 and the Maunder minima. 

We assume that the long-term changes of the solar irradiance 
are proportional to the changes of the solar activity 0(f) 



F^(t)-F Ll (t ) 0(f) -0(f o ) 



F LT (f M ) - F LT (f ) <p(t M ) - cf>(t y 



(A.4) 



where fo refers to our reference period (1976-1996 average) and 
fM refers to the Maunder minimum. 



Substituting F LT (f M ) - F LT (fo) = V LT F LT (t ), one can con- 
nect the relative changes of the solar flux and activity 



F LT (t) - F LT (f ) = 



y L TF LT (fo) 
<P(tm) - <f>(k)) 



■W)-<t>(to)). 



(A.5) 



Equation (|A3) links RMS 20 (F LT (f)) and RMS 2O (0(f)) via the pa- 
rameter Vlt- In combination with Eqs. ( lA.2b and dA.3b it allows 
one to reconstruct the solar variability backwards in time as a 
function of Vu and Vlt parameters. 

The approach, presented above is based on three assump- 
tions: 

1. The amplitude of the activity cycle is proportional to the 
long-term (i.e. cycle-averaged) solar activity. 

2. The long-term changes of the solar irradiance are propor- 
tional to the long-term changes of the solar activity. 

3. The long-term solar activity can be reconstructed from the 
cosmogenic isotope data. 

Let us note that these assumption are either employed or ful- 
filled because of the physical reasons in most of the available re- 
constructions, independently from the proposed mechanism for 
the secular changes. 

To estimate the accuracy of our approach we consider 
three different reconstruc tions of the solar irradiance to the 
past: lKrivovaetalJ (l2010l) (V LT =0.044%, Vu =0.000019 mag), 
Shapiro et al. (2011) with removed long -term trend (V r, T=0% , 
Vi i=0.0000 18 mag), and the original Shapi ro et all CO 1 11) 
reconstruction (V LT =0.38%, V n =0.000018 mag). They yield 
0.0001 mag, 0.00009 mag, and 0.00036 mag for the mean so- 
lar RMS 2 o variability over the 1620-1995 period. The same val- 
ues calculated with the set of Eqs. IA.214A.5l are 0.00013 mag, 
0.00012 mag, and 0.00035 mag. One can see that our approach 
slightly overestimates the first two numbers. The reason for this 
is that to allow a homogeneous reconstruction over the millennia 
we use the solar modulation potential to scale the cyclic vari- 
ability (see Eq. IA.3b instead of the sunspot number. The an- 
nual sunspot number is very close to zero during the Maunder 
minimum period, so the reconstructions predict the cessation of 
the cyclic variations (see Fig. |2j. At the same time the mod- 
ulation potential and consequently the cyclic variability in our 
approach did not go t o zero during the Maunder minimum 
(McCrack en et al.ll2004l) . Such deviations are comparable to the 
uncertainties of the reconstructions and are not essential for our 
conclusions. 

The calculation of the Ca II H and K lines, which are used 
to determine the solar chromospheric activity, made even more 
complicated by the effect of th e non-local thermodynamic al 
equilibrium (Ermol lTet al.l 12010). So the mean solar chromo- 
spheric activity cannot be directly deduced from the available 
reconstructions of the solar irradiance over the millennia . Instead 
we will linearly scale it with the long-term solar activity 



<P(t) 

0('o) 



((^Hk)o ~ (^HK)min) , 



(A.6) 



where (^^ K )o is the present value of t he mean solar chromo- 
spheric activity (log/?^ K =-4.895, (see lLockwood et all 120071) 
and (/^ K )mj n is the value of the mean chromospheric activity 
corresponding to = 0. 

In Fig. [3] we consider two extreme scenarios (see Sect. |3), 
one without the secular component in the solar chromospheric 
activity, and another with a strong secular component. Without 
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the secular component the changes of the mean solar chromo- 
spheric activity are only caused by the variable amplitude of 
the solar activity cycle and the (^ K )mj n value corresponds to 
the present value of the solar activity at minimum conditions 
(log^ K =-4.98; (see lJudge&Saarj ;2007). To describe the case 
of a strong secular component in the chromospheric activit y we 
put th e (R[j K ) m in value equals to the boundary value from Saar 
gOOi (log R' HK = -5.08). 

Appendix B: Possible constraints on the historical 
solar variability 

The amplitude of the secular solar photometric variability can 
be constrained by demanding the temporal mean solar variabil- 
ity be in agreement with the variability, indicated by the stellar 
data. The latter depends on the version of the variability versus 
activity regression and the scenario of the chromospheric activ- 
ity behaviour. To investigate different cases in Fig. IB.ll we plot 
the long-term solar variability Vlt as a function of the mean solar 
variability (RMS20) given by the stellar data, for two values of 
the 1 1 -year variabili ty V\ \ . The in clination effect was calculated 
according to lKnaack et al.l d2001l) . 

Let us note, however, that the positions of the Sun-like stars 
in Fig. [3] are also not necessarily fixed in time. This can affect 
the coefficients of the regression line and accordingly the esti- 
mated variability of the Sun. For example it is possible that, by 
coincidence, most of the stars among the twenty-one used for 
our analysis represent the periods of the relatively low (or high) 
variability. This will lead to the shifted position of the regres- 
sion line and consequently to the deviations in the V\ ( parameter 
determined using this limited selection of stars. 

If the behaviour of the stellar variabilities after the regres- 
sion to the solar level of the chromospheric activity is identical 
to the solar variability, i.e. the solar and regressed stellar variabil- 
ities can be parameterised employing the approach described in 
Sect. [3] and Appendix [A] then the resulting uncertainty can be 
found by performing the Monte Carlo simulation. To do this we 
considered a large number of sets consisting of 21 stars. Every 
star in these sets represents the Sun at some random twenty- 
year interval period (chosen from the full 9000-year dataset) and 
with a random angle between the stellar rotation axes and the 
direction to the observer. The RMS20 variability of these stars 
is calculated, emplo ying the approach described in Sect. [3] and 
AppendixlAland the lKnaack et alJ (1200 lh dependency of the 11- 
year variability on the angle between the solar rotation axes and 
the direction to the observer. For every set we calculated the 
mean RMS20 variability of the stars in the set, so for every pair of 
Vlt and V\\ parameters we have a distribution of mean RMS20 
variabilities. This allows us to calculate the range of Vlt and 
Vn parameters which can lead to the present value of the linear 
coefficients of the activity versus variability regression line. In 
Fig. IB. II we mark the resulting 2<x uncertainty in Vlt parameter 
with a shaded area. 

The mean solar log R' HK index in the case of the strong sec- 
ular component in chromospheric activity is -4.96 (see Fig. |3). 
Using the solid blue line in Fig. [T]one can find that this index 
corresponds to the 8.2 • 10~ 4 RMS variability. This value is in- 
serted in the left panel of Fig. IB. H and one can see that it yields 
Vlt = (0.43% - 0.76%). 

In the right panel of Fig. IB. II we also plotted the contour 
with Vn = 0.0002 mag and indicated four different values of the 
expected solar variability, corresponding to different treatments 
of the chromospheric secular changes and two regressions from 
Fig-Q] One can see that the minimum value of Vlt equals 0.17%. 



The parameter Vlt characterises the long-term variability as 
it is observed in the Stromgren b and y filters. Its connection with 
the variability at another spectral domain (or TSI) depends on the 
mechanism of the secular changes. For example a Vlt equal to 
0.17% leads to a 2.7 W/m 2 TSI change between the Maunder 
and th e last solar minima, adopting the model of Shapiro et al. 
(1201 lh . and to 1.9 W/m 2 , assuming that it is caused by the vari- 
ations of the solar effective temperature. 



A. I. Shapiro et al.: The place of the Sun among the Sun-like stars, Online Material p 4 




2 4 6 8 10 12 2 4 6 8 10 12 
RMS 20 , 10" 4 mag RMS 20 , 10" 4 mag 

Fig. B.l. Left panel: The dependency of the Vu parameter on the solar variability, retrieved from the stellar data. The shaded 
area indicates the uncertainty (2(f) due to the limited number of stars. The Vu parameter is set to 0.00044 mag. The dashed line 
corresponds to the expected mean solar variability RMS20 = 0.00082. The projection onto Vcr-axis yields the 95% interval for 
expected long-term variability of the Sun: 0.43% < Vu < 0.76%. Right panel: The same as the left panel, but the contour calculated 
with Vu = 0.0002 mag is added. The four dashed lines correspond to four values of expected mean solar variability RMS20: 
0.00097 mag (no chromospheric secular changes, regression with solid line from Fig. [TJ, 0.00082 mag (strong chromospheric 
secular changes, regression with solid line), 0.00065 mag (no chromospheric secular changes, regression with dashed line), 0.00053 
mag (strong chromospheric secular changes, regression with dashed line). The projection onto the Vcr-axis yields the interval 
0.17% < V LT < 1-2%. 
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Fig.B.2. Upper panel: RMS20 back to 400 BC, adopting V LT = 0.001 (lower curve) and V LT = 0.005 (upper curve). V u is set to 
0.0002 mag. Lower panel: Modulation potential back to 400 BC. The interval in the right upper corner denotes the satellite era of 
continuous solar irradiance measurements. 



